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We calculate the cross section of inclusive direct J/y/ photoproduction at next-to-leading order 
within the factorization formalism of nonrelativistic quantum chromodynamics, for the first time 
including the full relativistic corrections due to the intermediate 1 Sq\ , and 3 fj 8 ' color-octet 
states. A comparison of our results to recent HI data suggests that the color octet mechanism is 
indeed realized in J/y photoproduction, although the predictivity of our results still suffers from 
uncertainties in the color-octet long-distance matrix elements. 



X 



XVIII International Workshop on Deep-Inelastic Scattering and Related Subjects 

April 19-23, 2010 

Convitto delta Calza, Firenze, Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Direct J /iff photoproduction in at NLO in NRQCD 



Bernd A. Kniehl 



The factorization formalism of nonrelativistic quantum chromodynamics (NRQCD) [jl]] pro- 
vides a consistent theoretical framework for the description of heavy-quarkonium production and 
decay. This implies a separation of process-dependent short-distance coefficients, to be calculated 
perturbatively as expansions in the strong-coupling constant a s , from supposedly universal long- 
distance matrix elements (LDMEs), to be extracted from experiment. The relative importance of 
the latter can be estimated by means of velocity scaling rules; i.e. , the LDMEs are predicted to scale 
with a definite power of the heavy-quark (Q) velocity v in the limit v <C 1. In this way, the theoreti- 
cal predictions are organized as double expansions in a s and v. A crucial feature of this formalism 
is that it takes into account the complete structure of the QQ Fock space, which is spanned by the 
states n = 2S+l lJf with definite spin S, orbital angular momentum L, total angular momentum J, 
and color multiplicity a = 1,8. In particular, this formalism predicts the existence of color-octet 
(CO) processes in nature. This means that QQ pairs are produced at short distances in CO states 
and subsequently evolve into physical, color-singlet (CS) quarkonia by the nonperturbative emis- 
sion of soft gluons. In the limit v — > 0, the traditional CS model (CSM) is recovered in the case of 
5-wave quarkonia. 

Fifteen years after the introduction of the NRQCD factorization formalism [jl]], the existence 
of CO processes and the universality of the LDMEs are still at issue and far from proven, despite 
an impressive series of experimental and theoretical endeavors. The greatest success of NRQCD 
was that it was able to explain the J/\jf hadroproduction yield at the Fermilab Tevatron [0], while 
the CSM prediction lies orders of magnitudes below the data, even if the latter is evaluated at 
next-to-leading order (NLO) or beyond [Q]. Also in the case of J/\j/ photoproduction at DESY 
HERA, the CSM cross section significantly falls short of the data, as demonstrated by a recent 
NLO analysis [Q] using up-to-date input parameters and standard scale choices, leaving room for 
CO contributions [||]. Similarly, the J/yf yields measured in electroproduction at HERA and in 
two-photon collisions at CERN LEP2 were shown [§, 0] to favor the presence of CO processes. As 
for J/y polarization in hadroproduction, neither the leading-order (LO) NRQCD prediction 
nor the NLO CSM one [j3|] leads to an adequate description of the Tevaton data. The situation is 
quite similar for the polarization in photoproduction at HERA [Q]. 

In order to convincingly establish the CO mechanism and the LDME universality, it is an ur- 
gent task to complete the NLO description ofj/y hadro- ^ and photoproduction [Q, ^], regarding 
both J/y yield and polarization, by including the full CO contributions at NLO. While the NLO 
contributions due to the 1 Sq and 3 ,s[ 8 ' CO states may be obtained using standard techniques the 
NLO treatment of 3 pj 8 ' states in 2 — > 2 processes requires a more advanced technology, which has 
been lacking so far. In fact, the 3 pj 8 ' contributions represent the missing links in all those previous 
NLO analyses [g], g], and there is no reason at all to expect them to be insignificant. Specifically, 
their calculation is far more intricate because the application of the 3 pj 8 ^ projection operators to 
the short-distance scattering amplitudes produce particularly lengthy expressions involving com- 
plicated tensor loop integrals and exhibiting an entangled pattern of infrared singularities. This 
technical bottleneck, which has prevented essential progress in the global test of NRQCD factor- 



ization for the past fifteen years, was overcome for the first time in Ref. [10], which we review 
here. 

In direct photoproduction, a quasi-real photon y that is radiated off the incoming electron e 
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Figure 1: Sample diagrams contributing at LO (a and d) and to the virtual (b and e) and real (c and f) NLO 
corrections. 



interacts with a parton i stemming from the incoming proton p. Invoking the Weizsacker- Williams 
approximation and the factorization theorems of the QCD parton model and NRQCD [jl|], the in- 
clusive J/y photoproduction cross section is evaluated from 

da(ep^J/ V + X) = J £ f dxdyf r/e (x)f i/p (y)(^[n])da(Yi^cc[n}+X), (1) 

i,n ^ 

where fy/ e (x) is the photon flux function, fi/ p (y) are the parton distribution functions (PDFs) of 
the proton, (0 J /V[n]) are the LDMEs, and do{yi — > cc[n] + X) are the partonic cross sections. 
Working in the fixed-flavor-number scheme, i runs over the gluon g and the light quarks q = u,d,s 
and anti-quarks q. The Fock states contributing through the order of our calculation include n = 
3 5j 1 ', 1 Sq \ 3 S^\ 3 pfK Example Feynman diagrams for partonic LO subprocesses yi — > cc[n] +X as 
well as virtual- and real-correction diagrams are shown in Fig. [j]. 

We now describe our theoretical input and the kinematic conditions for our numerical analysis. 
We set m c = mj/y/2, adopt the values of mjiy, m e , and a from Ref. [11], and use the one-loop 

(two-loop) formula for aj . (ju), with «/ = 3 active quark flavors, at LO (NLO). As for the proton 
PDFs, we use set CTEQ6L1 (CTEQ6M) [0] at LO (NLO), which comes with an asymptotic scale 
parameter of aJJ d = 215 MeV (326 MeV), so that A^ D = 249 MeV (389 MeV). We evaluate 
the photon flux function using Eq. (5) of Ref. [pi] with the cut-off gLx = 2 Gey2 H3 on tne 



photon virtuality. Our default choices for the renormalization, factorization, and NRQCD scales 
are \i r = \if = m T and \i^ = m c , respectively, where m T = \j p\ + 4m^ is the J/y transverse mass. 
We adopt the LDMEs from Ref. Jl6|], which were fitted to Tevatron I data using the CTEQ4 PDFs, 
because, besides the usual LO set, they also comprise a higher-order-improved set determined by 
approximately taking into account dominant higher-order effects due to multiple-gluon radiation 
in inclusive J j\\f hadroproduction, which had been found to be substantial by a Monte Carlo study 
[g7|] . We disentangle ( G 1 ^ ( 1 S [ q ] ) ) and ( G ] l f ( 3 P^ ] ) ) , a linear combination of which is fixed by the 
fit only, as in Ref. [Q. The LO CO LDMEs are similar to the those obtained in Ref. Q by fitting 
Tevatron II data using the CTEQ6L1 PDFs [|l^]. The higher-order-improved CO LDMEs are likely 
to undershoot the genuine ones, which are presently unknown. 

Recently, the HI Collaboration presented preliminary data on inclusive J/ \j/ photoproduction 
taken in collisions of 27.6 GeV electrons or positrons on 920 GeV protons in the HERA II labora- 
tory frame [15]. They nicely agree with their previous measurement at HERA I [jilt]. These data 
come as singly differential cross sections in p\, W = a/ (Py + Pp) 2 , and z= (pj/y- Pp)/(py Pp), 
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(a) P?[QeV 2 ] (b) W [GeV] (c) z 



Figure 2: (a) p\, (b) W, and (c) z distributions of inclusive J/y/ photoproduction at LO and NLO in the CSM 
and full NRQCD in comparison with HI data [|4[ O]. The shaded (yellow) bands indicate the theoretical 
uncertainty due to the CO LDMEs. 



in each case with certain acceptance cuts on the other two variables. Here, p r , p p , and pjjy are 
the photon, proton, and // y four-momenta, respectively. In the comparisons below, we impose the 
same kinematic conditions on our theoretical predictions. 

The HI measurements [14, [l5|] of the p\, W, and z distributions of inclusive J /iff photopro- 
duction are compared with our new NLO predictions in full NRQCD in Fig. ^a)-(c), respectively. 
The uncertainty due the LDMEs is indicated by shaded (yellow) bands, whose upper margins (solid 
lines) refer to the LO set. For comparison, also the default predictions at LO (dashed lines) as well 
as those of the CSM at NLO (dot-dashed lines) and LO (dotted lines) are shown. Notice that 
the experimental data are contaminated by the feed-down from heavier charmonia, mainly due to 
y' — > J/\j/+X, which yields an estimated enhancement by about 15% [0]. Furthermore, our predic- 
tions do not include resolved photoproduction, which contributes appreciably only at 0.3 [|Trj|], 
and diffractive production, which is confined to the quasi-elastic domain at z ~ 1 and pj 0. These 
contributions are efficiently suppressed by the cut 0.3 < z < 0.9 in Figs. §(a) and (b), so that our 
comparisons are indeed meaningful. We observe that the NLO corrections enhance the NRQCD 
cross section, by up to 115%, in the kinematic range considered, except for z<;0.45, where they 
are negative. As may be seen from Fig. |2](c), the familiar growth of the LO NRQCD prediction 
in the upper endpoint region, leading to a breakdown at z = 1, is further enhanced at NLO. The 
solution to this problem clearly lies beyond the fixed-order treatment and may be found in soft 
collinear effective theory [20]. The experimental data are nicely gathered in the central region of 
the error bands, except for the two low-z points in Fig. ||(c), which overshoot the NLO NRQCD 
prediction. However, this apparent disagreement is expected to fade away once the NLO-corrected 
NRQCD contribution due to resolved photoproduction is included. In fact, the above considera- 
tions concerning the large size of the NLO corrections to hadroproduction directly carry over to 
resolved photoproduction, which proceeds through the same partonic subprocesses. On the other 
hand, the default CSM predictions significantly undershoot the experimental data, by typically a 
factor of 4, which has already been observed in Ref. [Q]. Except for p?^4 GeV 2 , the situation is 
even deteriorated by the inclusion of the NLO corrections. 

Despite the caveat concerning our limited knowledge of the CO LDMEs at NLO, we con- 
clude that the HI data [ 14, 15] show clear evidence of the existence of CO processes in nature, as 
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predicted by NRQCD, supporting the conclusions previously reached for hadroproduction at the 
Tevatron [§] and two-photon collisions at LEP2 [Q]. In order to further substantiate this argument, 
it is indispensable to complete the NLO analysis of inclusive J/y hadroproduction in NRQCD, 
by treating also the 3 pj 8 ' channels at NLO, so as to permit a genuine NLO fit of the relevant CO 
LDMEs to Tevatron and CERN LHC data. This goal is greatly facilitated by the technical advance- 
ment achieved in the present analysis. 
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